Stochastic growth of quantum fluctuations during slow-roll inflation 
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We compute the growth of the mean square of quantum fluctuations of test fields with small effec- 
tive mass during a slowly changing, nearly de Sitter stage which takes place in different inflationary 
models. We consider a minimally coupled scalar with a small mass, a modulus with an effective 
mass oc H 2 (with H the Hubble parameter) and a massless non-minimally coupled scalar in the test 
field approximation and compare the growth of their relative mean square with the one of gauge 
invariant inflaton fluctuations. We find that in most of the single fields inflationary models the mean 
square gauge invariant inflaton fluctuation grows faster than any test field with a non-negative ef- 
fective mass. Hybrid inflationary models can be an exception: the mean square of a test field can 
dominate over the gauge invariant inflaton fluctuation one on suitably chosen parameters. We also 
compute the stochastic growth of quantum fluctuations of a second field, relaxing the assumption 
of its zero homogeneous value, in a generic inflationary model; as a main result, we obtain that the 
equation of motion of a gauge invariant variable associated, order by order, with a generic quantum 
scalar fluctuation during inflation can be obtained only if we use the number of e-folds as the time 
variable in the corresponding Langevin and Fokker-Planck equations for the stochastic approach. 
We employ this approach to derive some bounds for the case of a model with two massive fields. 

PACS numbers: 04.62,+v, 98.80.Cq 



I. INTRODUCTION 

The theory of quantum fields in an expanding universe 
has evolved from its pioneering years [l| into a necessary 
tool in order to describe the Universe on large scales. The 
de Sitter background - characterized by the Hubble pa- 
rameter H = a/a being constant in time (for the flat spa- 
tial slice), where a(t) is the scale factor of a Friedmann- 
Robertson- Walker (FRW) cosmological model - has been 
the main arena to compute quantum effects even before 
becoming a pillar of our understanding of the early infla- 
tionary stage and of the recent acceleration of the Uni- 
verse. 

However, while \H\ <C H 2 for any inflationary model, 
H may not become zero in a viable model, apart from 
some isolated moments of time. Indeed, the standard 
slow-roll expression for the power spectrum of the adia- 
batic mode of primordial scalar (density) perturbations 
becomes infinite, i.e., meaningless, if H becomes zero dur- 
ing inflation.* Outside the slow- roll approximation, H 
may reach zero [3|, but for a moment only. Therefore, 



This is why the statement sometimes found in literature that, 
for H = const, a flat (Harrison-Zeldovich) n s = 1 spectrum 
is generated is also meaningless. Actually, it is a V(<f>) oc </>~ 2 
inflaton potential that leads to n B = 1 (and r 7^ 0) in the slow- 
roll approximation, see for the exact solution for V(rp) without 
using this approximation. 



the study of quantum effects in a nearly de Sitter stage 
with H 7^ 0, in particular, when the total change in H 
during inflation is not small compared with its value dur- 
ing the last e-folds of inflation 0-0] (see also the recent 
papers [§-flol]). is not of just pure theoretical interest. 
Among the main results of previous investigations, it has 
been shown that the infrared growth of minimally cou- 
pled scalar fields with a zero or small mass m <C H in a 
background with a practically constant H [TT| - tl3| occurs 
for massive fields when H changes significantly during 
inflation [H, @, and that the stochastic approach, orig- 
inally mainly applied to a new inflationary type back- 
ground with a small change in H during inflation flij . 
also works in realistic chaotic type inflationary space- 
times 0] (here we do not discuss its application to eternal 
inflation 3.3 and to interactions in a de Sitter back- 
ground [imi) . a particle production in a realistic infla- 
tionary background is so different from the corresponding 
one in the de Sitter space-time that it has prompted us 
to reconsider the amplification of nearly massless mini- 
mally coupled scalar fields in inflation with a quadratic 
potential 0, @] and compare it with the dynamics driven 
by a scalar field condensate. 

In this paper we wish to tackle in more detail the mod- 
uli problem issue. On one hand we wish to extend our 
results in [7] to different inflationary models and to dif- 
ferent types of test fields, not only to massive minimally 
coupled scalar fields, but also to massless non-minimally 
coupled scalar fields and moduli with an effective mass 
oc H 2 . We investigate effects 0(H) and we therefore need 



2 



to consider both these extensions. On the other hand, it 
is known that the mean square of gauge invariant infla- 
ton fluctuations grows [Ef and the stochastic description 
for this effect for a general potential has been established 
7]; it is therefore interesting to compare the quantum 
amplification of test fields not only with the background 
inflaton dynamics, but also with the stochastic growth 
of gauge invariant inflaton fluctuations. This compari- 
son aims for a self-consistent understanding of quantum 
foam during inflation. 

We then discuss the diffusion equation for general 
scalar fluctuations in a generic model of inflation. On 
using the results obtained by field theory methods, we 
show that the diffusion equation for the gauge invariant 
variable associated with this generic scalar fluctuation 
should be formulated in terms of the number of e-folds 
N. 

The paper is organized as follows. In Section II we 
choose four representative cases of the inflationary "zoo " 
on which we focus in this paper. In Section III we 
compute the stochastic growth of a minimally coupled 
scalar with a small mass, a modulus with an effective 
mass cx H 2 , and a massless non-minimally coupled scalar 
in the test field approximation for the four inflationary 
models considered. In Section IV we review the result 
obtained in Ref. Q for gauge invariant inflaton fluctua- 
tions and compare it with the growth of the test fields in 
the four different inflationary models considered. In Sec- 
tion V we examine the stochastic approach for a two field 
model for two generic self-interacting potentials, and we 
compare our solution, obtaining some constraints on the 
parameters, for a particular two quadratic field model. 
In Section VI we illustrate our conclusions. 

Our paper does not include a derivation of the stochas- 
tic approach. Indeed an introduction of the stochastic 
method is given in Refs. [14,18] and a comparison be- 
tween stochastic methods and quantum field theory re- 
sults is done in Refs. [14,18] and in our previous paper 
[7]- 



II. INFLATIONARY MODELS 



During the slow-roll trajectory we have: 
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H? - -m 2 N 



mM, 



P ] ■ 



(2) 

(3) 
(4) 



where M nl = 8nG (these formulas were obtained al- 
ready in [2(| in the context of a closed bouncing FRW 
universe with two quasi-de Sitter stages during contrac- 
tion and expansion). Let us note that in the numerical 
results presented in the figures all dimensional quantities 
have been rescaled w.r.t. m p i = Mpi^/Sn. We then con- 
sider the case of a quadratic potential (of arbitrary sign) 
uplifted with an offset Vq: 



v(4>) 



(5) 



With the positive sign the potential in Eq. ([5]) is an 
approximation for the simplest model of hybrid inflation 
well above the scale of the end of inflation; in this case 
<f> decreases during the inflationary expansion. With the 
negative sign the potential in Eq. ([5]) is a simple small 
field inflation model, again far from the end of inflation; 
in this case <fi increases during the inflationary expansion. 
In the following we use the approximate solution for the 
square of the inflaton as: 



2 e T2- 



- (N—Ni) 



(6) 

The exact expression can be given in the implicit form 



valid when H 



V Q /(UP pl ), i.e. M_|^ 2 



Vo 



2M 2 M 



IM 2 



N-N t = T 

'"pi ~~ P i 
As another large field inflationary model we consider 
an exponential potential 



The detailed evolution of the expansion during the ac- 
celerated stage depends on the inflaton potential and so 
does the growth of quantum fluctuations. For this rea- 
son we consider in the following four different potentials 
which are representative of the "inflationary zoo ". Since 
we shall study the growth of quantum fluctuations as a 
function of the number of e-folds 



N = log 



a(U) 



(1) 



V = V e . (7) 

This potential leads to a power-law expansion given by: 



(8) 



with p = 2/A 2 [2l| (we consider a(ti) = 1 for the sake 
of simplicity). Such a solution is stable for p > 1 and 



we shall give the evolution of the Hubble parameter as a 
function of N. 

The first obvious model is chaotic quadratic inflation, 
which we have also used in our previous investigations 



t For a double well SSB potential V = \(<f> 2 — <f> 2 ) 2 , which may- 
be approximated by Eq. J5]l for small field values, the condition 
becomes \(<j> 2 - <j> 2 )/4> 2 \ < 1 together with (f> 2 /4> 2 < 1, which 
again means that V{4>) is always considered far from the minima. 
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the slow-roll conditions are well satisfied for p >> 1. In 
particular one obtains 



H = Hi e- N 'v = * e~ N /P 
U 



M p iN + fa 



(9) 
(10) 



which is just the integral form of Eq. (13) of Ref. Q 
generalized to an arbitrary inflaton potential. 

For the quadratic inflaton case we report here the so- 
lution given in Ref. Q § : 



(x 2 )ren 



3ir 



87r 2 (2m 2 — mt) 



4-2- 



-H 



4-2- 



(13) 



III. GROWTH OF TEST FIELDS WITH SMALL 
EFFECTIVE MASS IN THE STOCHASTIC 
APPROACH 

We shall consider a test scalar field with a small effec- 
tive mass and a zero homogeneous expectation value on 
an inflationary background driven by an inflaton with po- 
tential V{fa) in the slow- roll approximation. The evolu- 
tion equation for the renormalized mean square (x 2 )ren 
(the pedix REN will denote renormalized in the follow- 
ing) in the next three subsections (Eqs. 11, 16 and 
22) follows in a straightforward manner from our pre- 
vious paper [7]. Note that the right hand side of Eqs. 
(11,16,22) representing the contribution of created fluc- 
tuations ("particles") is obtained under the natural as- 
sumption of the absence of particles in the in-vacuum 
state, more exactly that each Fourier mode k of the quan- 
tum field x w & s hi the adiabatic vacuum state deep inside 
the Hubble radius and long before the first Hubble radius 
crossing during inflation, i.e., when its energy uj = k/a(t) 
was much larger than H(t). * So, the explicit time- 
asymmetry of Eq. (11) shows that this in-vacuum is not 
de Sitter invariant; it is unstable and creation of fluctu- 
ations (particles) of light scalar fields, as well as metric 
perturbations, takes place. In turn, the cause of this in- 
stability may be finally traced to the expansion of the 
Universe. 



where we have assumed (x 2 )REN(^i) = (we shall adopt 
the same choice afterwards if not otherwise stated). We 
then consider the potentials in Eq. ([5]), in the lowest 
order approximation; that is, for V ~ V = SHqM 2 ^ we 
have 



(x 2 



'REN 



3^ 
8ir 2 mi 



1 



-JV 



(14) 



Let us note that the corrections induced by a non-zero 
M 2 4> 2 /Vo term are typically small both for the case of 
hybrid inflation as well as for small field inflation (as 
long as the field does not grow too much due to instabil- 
ity). The corresponding analytic expressions, obtained 
using Eq. ©, can be written in terms of hypergeomet- 
ric functions but we do not report them here. For the 
exponential potential we obtain: 



(x 2 )ren 



P 2 

■^Hi exp 





3 Hf 
■-—%Ei 



3H 2 
3H 2 



exp Uw 



(15) 



where E% is the exponential integral function (see, for 
example, [Hj|). 



Growth of scalar fields with m x in the 
stochastic approach 



The stochastic equation is: 



d(x 2 >REN , 2 ™ 2 



dN 3H 2 (N) 
Its general solution is 



(x 2 )ren = (x 2 )ren(^) 



-dn 



x (x 2 )ren 



H 2 (N) 
An 2 



N 



dn t ' e J 
An 2 



3H 2 {fi 



(11) 



■d,7i 



(12) 



B. Growth of moduli fields with 
stochastic approach 



cH in the 



If |c| <C 1, the stochastic equation takes the form: 
d(X 2 >REN , 2c. 2N H 2 (N) 



dN 3 
Its general solution is 



(x 2 )ren 



Air 2 



(x 2 



'REN 



(X 2 )REN(iV l ) + 



N 



dn \ 
An 2 



(16) 



e 3 



(17) 



: 2m 2 we obtain 

+ This assumption also means the absence of the so called trans- 
Planckian particle creation, see |2'J and references therein for (x 2 )ren — ] Q g ( Ml. 

discussion why trans-Planckian particle creation should be ab- 4ir 2 m 2 \ H 

sent in the standard QFT. 
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For V(4>) = m <fr/2 we obtain: 



(x 2 ! 



REN 



m 

6V 2 

3H? 
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4? 



2c 



-Nn 



2c 



-N 



(18) 

where Nt = is equal to maximal number of possible 
e-folds in this chaotic model. In the limiting case c — > 
and at the end of inflation (JV = Nt — 3/2), we recover 
the result 1 : 



(x 2 



'REN 



12tt 2 



-Nl = 



3Hf 



I6n 2 m 2 



(19) 



For the potential in Eq. ([5]), we consider the lowest 
approximation as in the previous subsection. Therefore 
the result can be simply obtained by substituting m 2 = 
cH 2 in Eq. (fill): 



(x 2 )ren — 



87T 2 CV 



fcTV 



For the power-law inflation case we obtain: 



(x 2 



'REN 



8tt 2 



(20) 



(21) 



Let us note that m 2 — V<f,<f, = |(1 — ^)H 2 for this 
particular model, so the results above are also valid for 
the case m 2 = cm 2 with c — §(1 — ^-) _1 c. 

™ This corresponds to the massless limit of moduli production com- 
puted in Eq. (15) of 7] for a -> 0. 



Growth of non-minimally coupled scalar fields 
in the stochastic approach 



The stochastic equation is now: 



d(x 2 



'REN 



dN 



+ 4£(2-e)( X 2 >REN = 



H 2 {N) 
4tt 2 



(22) 



where £ is the non-minimal coupling to the Ricci scalar R 
and we assume that |£| <C 1 (however, £iV may be large). 
Indeed, the term in the action proportional to £x 2 -R gives 



an effective time dependent mass for \: 
where e = —jj2- 



6efP(2-e) 



Its general solution is 



(x 2 )ren — 



(x 2 )ren(^) 



JY 



17-2+45/- \ 

dn fl — w e « n 

4n 2 H?Z 



H(N)J 



-S£N 



(23) 



If we again consider the chaotic scenario induced by a 
massive inflaton field as in the previous subsections, the 
integral can be easily computed in a closed form in terms 
of the exponential integral function E v (z). Assuming 
initially, one finds 



(x 2 )ren 



(x 2 )ren 



m 2 e 8ttN T ~N) r 



6tt 2 (Nt-N) 



(Nt~N) 2 + 2 ^E^ 2 ^(N t ~N)) - (NT-Nf+^E^^faiNT-N)) 



(24) 



One can verify that in the limit £ —> at the end of 
inflation and a fixed large value for Nt, the result of 
Eq. (fll?)) for a massless modulus is again reobtained. 

For the potential in Eq. ([3]) we have, in the same ap- 
proximation as in the two previous subsections, 



(x 2 )ren 



32tt 2 ^ V 



-8£iV 



(25) 



As before, let us finish with the case of a power-law 
model of inflation. We obtain: 



(x 2 )ren 



p 

8tt 2 



Hf (-2£ - 1 



4Kr' (e 



(26) 



IV. COMPARISON WITH THE GROWTH OF 
INFLATON FLUCTUATIONS 

The results of the previous section should be com- 
pared with the growth of gauge-invariant inflaton fluc- 
tuations 5(f>, the Mukhanov variable [24j which is used 
to canonically quantize the Einstein-Klein-Gordon La- 
rangian. The evolution equation for (<5</> 2 )ren found in 
can be re-written as: 



d(<5</> 2 )REN 



dN 



+ 2 (77 - 2e) 



;ren 



H 2 (t) 
At: 2 



(27) 



where 



M 2 J-f 



(28) 
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In Eq. (|27p the positivity of rj — 2e is not determined 
by the convexity of the potential, i.e. > 0, as we 

would expect in the absence of metric perturbations. The 
threshold corresponds to the following condition on the 
potential: 



> 0. 



(29) 



With the use of the slow-roll expressions for the scalar 
spectral index n s and the tensor-to-scalar ratio r, one 
obtains: 
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n s — 1 = — 6e + 2?y 
r = 16e, 

and Eq. (|27fl can be rewritten as: 



(30) 
(31) 



d(S(/) 2 



•hen 



dN 



(n s -l + (S4> 2 ) 



REN 



H 2 (t) 
4tt 2 



Both Eqs. (|29|) and (|32|) tell us that power-law inflation, 
for which n s — 1 = — r/8 holds, lies at the threshold be- 
tween two opposite behaviors. Power-law inflation with 
78 < p < 246 is allowed at the 95 % confidence level [H[ . 
We note that Eq. (|32|) is the same for a modulus with 
the mass m 2 = cH 2 and c = 3(n s — 1 + r/8)/2: below 
the power-law inflation line inflaton fluctuations behave 
as a modulus with negative c. 




FIG. 1: Evolution of the mean square quantum fluctuations 
(in units of mp^) versus the number of e-folds N for the 
quadratic chaotic model. For the inflationary background 
we have chosen the inflationary trajectory in Eq. <[3j) with 
m = 10 -6 m p i and Hi = 10 m. The mean square gauge in- 
variant inflaton fluctuation (thick line) dominates over those 
of test fields (m x ~ 0.3m is the solid line, c = 0.02 is the 
dashed line, £ = 0.001 is the dotted line). 



2: Evolution of the mean square quantum fluctuations 
units of m pl ) versus the number of e-folds N for the 
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small field inflationary model in Eq. (0). 
tionary background we have chosen Vo = 



For the infla- 



2.6 x 10" 



'pi- 



M — 0.85 x 10 m p i and 4>i — 0.3 m p i as parameters. The 
mean square gauge invariant inflaton fluctuation (thick line) 
(32) dominates over those of test fields (m x = W~ 2 Hq is the solid 



line, c = 0.1 is the dashed line, £ = 0.05 is the dotted line). 



For the potential in Eq. (|5|) we obtain the following 
in the lowest non-trivial approximation for small in- 
dependent corrections in the potential: 



/REN 



± 



V' 2 
1 o 



24tt 2 M 2 M4 



Vo 



-(N-Ni) 



(35) 

where, as discussed previously, the case with a minus sign 
in the exponent refers to the hybrid model whereas the 
other case is associated with the small field inflationary 
model. Let us note that for a small value of the ex- 
ponent (an almost constant potential) or for very small 
(N — iVi), by expanding up to the linear order, one ob- 
tains almost the case of a dc Sitter background (Al 2 = 0), 
with (50 2 )ren linearly growing in N. In this approxima- 
tion we see that the hybrid model is characterized by 
(<50 2 )ren bounded as N — > oo. This statement remains 
true even after dropping our approximations (see below) . 

This leads, by invoking the consistency of the pertur- 
bative expansion in field fluctuations through the condi- 
tion (<5</> 2 )ren "C 4> 2 , to the following hierarchy which the 
inflationary model has to satisfy: 



The solution of Eq. (J17J) is: 



4tt 2 



, 2\ <N) r H 2 (n) 

b )ren = -t— r- / an- 



e(n) 



(33) 



Vo 



247r 2 M 4 ! 



M 2 cj) 2 1 

<< ^T <<1 



(36) 



For the quadratic chaotic potential the solution was 
found in 



/REN 



Hf - H 6 
8n 2 m 2 H 2 



(34) 



We have also computed the expression for the fluctua- 
tions by solving Eq. (|2T|) using the expression in Eq. © 
with no further approximations. In this more general 
case we obtain 
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4F 2 (1 - y) + 3M*<ft y UM^N N t ) + 0f(l - y)) ± 2/(1 - y 2 )^ 
i 2 )REN-± ^ Af2 ^ 2 ' (37) 



2M 2 

where we have set y — y(N) — e T ^ ^ N N ^ . From 
this expression, when analyzing the hybrid inflation case, 
one can notice that the fluctuations have a maximum for 
a certain amount of e-folds N and then decay to the 
asymptotic value for a large number of c-folds. Never- 
theless, such a maximum is typically a few percent above 
the asymptotic value which has been already obtained 
above using a more crude approximation. 
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FIG. 3: Evolution of the mean square quantum fluctuations 
(in units of m 2 ^) versus the number of e-folds N for the hy- 
brid model in Eq. ([5]). For the inflationary background we 
have chosen Vo = 2.6 x 10" 12 m* b M — 1.8 x 10~ 6 m p i and 
4>i = 0.3 m p i as parameters. In this case the mean square of 
moduli can dominate over the mean square of gauge invariant 
inflaton fluctuation (thick line): the parameters chosen are 
m x = 10~ 2 iio (solid line), c = 0.002 (dashed line), £ = 0.05 
(dotted line). 

For the exponential potential we find 

W 2 )REN-^(^-ff 2 ), (38) 

which at late times (see also (26[) becomes: 

(^ 2 )ren = |J • (39) 

We show in Figures [IJ [2 [3] and HI one for each infla- 
ton potential investigated, the mean square of quantum 
fluctuations of the three types of test fields together with 
the gauge invariant inflaton ones. For the hybrid model, 
quantum fluctuations of test fields with a small effective 
mass can dominate the gauge-invariant inflaton one, be- 
cause of the presence of the leading constant term in the 
potential. 
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FIG. 4: Evolution of the mean square quantum fluctuations 
(in units of TUpi) versus the number of e-folds N for the ex- 
ponential potential. For the inflationary background we have 
chosen the inflationary trajectory in Eq. ((9| with p — 100 
and ti = 10 wipj 1 . The mean square gauge invariant infla- 
ton fluctuation (thick line) dominates over those of test fields 
(m x = 10 -6 m p i is the solid line, c = 0.1 is the dashed line, 
£ = 0.05 is the dotted line). 



V. GROWTH OF QUANTUM FLUCTUATION 
IN TWO FIELD INFLATIONARY MODELS 

We now wish to consider a two field model in which 
an inflaton </> and a minimally coupled scalar field x are 
present (see [13] for a different approach to the mod- 
uli problem). We shall neglect the \ ener gy density 
and pressure in the background FRW equations. We 
expand to second order in the uniform curvature gauge 
(UCG), in which the inflaton fluctuation ip coincides with 
the gauge-invariant Mukhanov variable, the Einstein and 
Klein-Gordon equations. 

In the test field expansion x{x, t) = Xo(t) + (%, t) + 
X^ 2 \x,t) + the homogeneous term satisfies 



Xo + 3H X0 + V x = , (40) 

while fluctuations satisfy, order by order, for the lead- 
ing order in the slow-roll approximation and in the long- 
wavelength limit (neglecting vector and tensor contribu- 
tions), the following equations: 



3H X ^+V xx x^=2^V x ^\ (41) 
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3H X {2) + V xxX {2) 



'XX/ 

2§ V x ^ 



+ 



H 



Vx^ (1)2 



XXX (1)2 

2 X 



(42) 



Following the consideration in section VI of Ref . , we 
wish to investigate which time variable in the stochastic 
equation should be chosen to re-obtain, order by order, 
the equation of motion for the test field x starting from 



dx 
dt 



V x 



(43) 



and expanding order by order. For a general time vari- 
able t = J H(t) n dt, the equation becomes 



dx 



H{t) n dt m{(t>) n+1 



(44) 



As before, expanding to leading order in the slow-roll 
approximation, we obtain the following equation to the 
first and second order: 



dx ( 



i) 



dt 



-—V, v* 1 ' 4- -(n 
3H v xxX + 3 {n + L) H2 V x ip 



(45) 



d X {2) 
dt 



^*x (2) +^+ 1 )fl^ (2) - 



6H 



xx V {1) X {1) — Q {n+l)V x 



r 



H, 



H 2 



+ (n 



H 3 
_J46} 



(1)2 



As is easy to verify, we recover the former result only for 
n = 1. So for the case of a test scalar field, evolving in 
a FRW inflaton driven space-time, in the UCG the right 
time to consider is the number of e-folds N = J H(t)dt, 
this recovers the result obtained in [7j for the inflaton 
fluctuations. As for the case of the standard Mukhanov 
variable Q, which is defined, order by order, as the value 
of the inflaton perturbation in the UCG, we can define 
a generic gauge-invariant Mukhanov variable Q x , associ- 
ated with the perturbation of Xi as the value that this 
perturbation has in the UCG. In this way Q x — x^ m 
the UCG and, as for the variable Q in [7J, the equations 
above can be regarded as the gauge- invariant equations of 
motion, to first and second order, of this new Mukhanov 
variable, where one replaces x*™' with and tp^ with 
Q (n) . 

As for the case of the Mukhanov variable, this result 
can be considered as a starting point to study the fluc- 
tuations of x in the stochastic approach for an arbitrary 
potential V in the described background. The correct 
stochastic differential equation is obtained with respect 
to the number of e-folds (dN = H(t)dt) which appears 
to be the right evolution parameter. One starts from 



the slow-roll approximation to the Heisenberg equation, 
which can be interpreted in a general non-perturbative 
sense, for the large-scale quantum field x 



dN^ 

(/ x (A 1 ,x 1 )/ x (7V 2 ,x 2 )) 



f,S(N^N 2 f n ^- X f 
4tt 2 xi -x 2 



where f x is the stochastic noise term given, to the leading 
order in the slow-roll approximation, by 



/ x (t,x) = eaH 2 



d 3 k 
(2^]3A 



■S(k - eaH) b k Xk{t)e 



+b{xt(t)> 



-Hk-x 



(47) 



Thus, on expanding to the second order, one obtains the 
following stochastic equations for x and x^ '■ 



dx^ 

dt 



— V Y (1) 



3H 



^ff^ (1) +/x, (48) 



d X {2) 
dt 



1 

3H 



V v (2) 



2 H ^ 
3IP 



6H 



v xxxX 



(1)2 



-Ely w (1) y (1) --V 
3H 2 V *xV X 3 V X 



Hi 



(1)2 



H, 



H 



Jx 



if- 



Let us consider the first order stochastic equation. Its general solution with the zero initial condition is given 



8 



by Taking into account that {ip^ 1 f x ) — 0, it is easy to de- 

rive the expression for the mean square of the first order 

A) =vj t l + k)dr. (50) fluctuation ((*X (1) ) 2 ) ■ 



r 



3H 2 ^ ' V x 



(X 



(1)2\ 



V£ J dT J d V 



^L dT 



9 H{t) 2 H{-qY 
H(t) 3 4 



V x (t) V x (tj) 



(fx(r)fx(v)) 



(51) 



where we have used the stochastic solution 



I* 
V 



r 



(52) 



with being the stochastic noise for the inflaton defined analogously to f x . Similarly, one can obtain the following 
solution for a vacuum expectation value of the second order fluctuation: 



(x (2) > = v x 



dT 



2^/ (2)) _ _LYxxx , (i )2) lE±Ymi^) x w\ + 1 ( Ei± _ 3 ^ /^a^ 



(53) 



where, to expand further, we should substitute Eq. (l5"Tj) . the result for ((p^) and (y' 1 ^ 2 ) obtained in Q and 



<^ (1) x (1) ) 



12lT 2 ffM^ Jt 



t pt 
dr I dr/ 



H(t) 5 H(rj) 

h(t) H( v r 



(54) 



r 



A. A working example: two field quadratic model 

2 j 2 

Let us now consider the particular case V(cf>) = — # 

- m 2 2 

and V(x) = — | — ■ Classical slow-roll inflation in this 
model and the evolution of small perturbations in it were 
calculated in [28| , but here we take into account the back- 
reaction of the generated quantum fluctuations of these 
scalar fields on the evolution of their background values. 
By solving the background equations, one obtains the 
following zero order solution for the test field x 



X < 0) (i)=X W (*i) 



(0), 



H(t) 
H{U) 



(55) 



It remains a test field for the whole duration of the infla- 
tion era if 



x i0 Hu) 2 « 



a m 

1 + 5F 



for any value of H (where a = — §■ ) . For the case a«l, 
we obtain the following limiting condition at the end of 
inflation (H ~ m): 



6 



a. 



(57) 



and for this particular background we can solve Eq.([5T 
obtaining 



(X 



(1)2\ _ 



3H 2 



8ir 2 m 2 (2-a 



H 4 



48tt 2 



H\ 1 , 2 

a TP {H 



ffif -(58) 



Thus, one obtains the term already considered in Q plus 
a new term which depends on the background value of \. 
At the end of inflation, the leading value of this second 
term is negligible with respect to the leading value of the 
first one, for a < 2, if 



x iQ Hu) 2 « 



18 1 M 2 im 2 
7 ^~d 2 H? 



(59) 



This condition is different from, and can be stronger 
than, the condition (|57j) . If we consider the particular 
case a< 1 and require that ([57]) implies (I5T)1) , we obtain 
the following condition on a: 



a < 



2H 2 



(60) 
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Analogously we can evaluate Eq . ([53 



(X 



(2)\ 



8tt 2 Af?, 



H 



H 

2 1 + a 

12 



#? 1 - a/2 Hfl-a 



H 4 6 



H 2 4 



(61) 



which to leading order gives 



(2)v g X (0) fe) 



48?r 2 Mi 



Hf \1 — 



(62) 



VI. CONCLUSIONS 

Motivated by previously found differences for gravita- 
tional particle production in the de Sitter background 
and in realistic inflationary models with H ^ 0, we have 
studied in detail the growth of quantum fluctuations for 
the latter case. We have selected four different poten- 
tials as representative examples of the inflationary zoo 
and different types of nearly massless fluctuations, in- 
cluding inflaton ones. We have rewritten in Eq. (|32[) the 
diffusion equation for the gauge invariant inflaton fluc- 
tuations found in Q, emphasizing the role of the slope 
of the spectrum of curvature perturbations n s and of the 
tensor-to-scalar ratio r, i.e. the relevant observable quan- 
tities. 

We have found that for most of the inflationary models, 
the mean square of the gauge invariant inflaton fluctua- 



tions dominates over those of moduli with a non-negative 
effective mass. Hybrid inflationary models can be an ex- 
ception: the mean square of a test field can dominate 
over that of the gauge invariant inflaton fluctuations on 
choosing parameters appropriately. Our findings show 
that the understanding of inflaton dynamics including its 
quantum fluctuations is more important than the moduli 
problem in most of the inflationary models. 

We have then discussed the stochastic approach for 
general scalar fluctuations, which may have a non zero 
homogeneous mode, in a generic model of inflation. We 
show, by using the field theory results as a guideline, that 
the stochastic equations for the gauge invariant variable 
associated with such scalar fluctuations are naturally for- 
mulated as a flow in terms of the number of e-folds N. 
Finally we have studied the particular case of a massive 
inflaton and a second massive scalar field Xj for which we 
show how to extract some bounds for the homogeneous 
mode of x an d for its mass. 



Acknowledgements 

G.M. was supported by the GIS "Physique des Deux 
Infinis". A.S. was partially supported by the Russian 
Foundation for Basic Research, Grant No. 09-02-12417- 
ofi-m. This work was started during his visit to Bologna 
in 2009 financed by INFN: we thank INFN for support. 



[1] N. D. Birrell and P. C. W. Davies, Quantum Fields in 
Curved Space (Cambridge University Press, Cambridge, 
1982). 

[2] A. A. Starobinsky, JETP Lett. 82, 169 (2005). 
[3] O. Seto, J. Yokoyama and H. Kodama, Phys. Rev. D 61, 
103504 (2000). 

[4] F. Finelli, G. Marozzi, G. P. Vacca and G. Venturi, Phys. 

Rev. D 65, 103521 (2002). 
[5] F. Finelli, G. Marozzi, G. P. Vacca, and G. Venturi, Phys. 

Rev. D 69, 123508 (2004). 
[6] F. Finelli, G. Marozzi, G. P. Vacca and G. Venturi, Phys. 

Rev. D 74, 083522 (2006). 
[7] F. Finelli, G. Marozzi, A. A. Starobinsky, G. P. Vacca 

and G. Venturi, Phys. Rev. D 79, 044007 (2009). 
[8] Y. Urakawa and K.-I. Maeda, Phys. Rev. D 77, 024013 

(2008); Phys. Rev. D 78, 064004 (2008). 
[9] T. M. Janssen, S. P. Miao, T. Prokopec and 
R. P. Woodard, Class. Quant. G rav. 25, 245013 ( 2008). 
[10] T. M. Janssen and T. Prokopec. larXiv: 0906. 06661 
[11] A. D. Linde, Phys. Lett. B 116, 335 (1982). 
[12] A. A. Starobinsky, Phys. Lett. B 117, 175 (1982). 
[13] A. Vilenkin and L. Ford, Phys. Rev. D 26, 1231 (1982). 
[14] A. A. Starobinsky, in Field Theory, Quantum Gravity 
and Strings, eds. H. J. De Vega and N. Sanchez, Lect. 
Notes in Physics 246, 107 (Springer, New York, 1986). 
[15] A. D. Linde, Phys. Lett. B 175, 395 (1986). 



[16] A. D. Linde, D. A. Linde and A. Mezhlumian, Phys. Rev. 

D 49, 1783 (1994). 
[17] N. C. Tsamis and R. P. Woodard, Nucl. Phys. B 724 

(2005) 295. 

[18] A. A. Starobinsky and J. Yokoyama, Phys. Rev. D 50, 
6357 (1994). 

[19] E. O. Kahya and V. K. Onemli, Phys. Rev. D 76, 043512 
(2007). 

[20] A. A. Starobinsky, Sov. Astron. Lett. 4, 82 (1978). 
[21] F. Lucchin and S. Matarrese, Phys. Rev. D 32, 1316 
(1985). 

[22] E. W. Kolb, A. A. Starobinsky and I. I. Tkachev, JCAP 

0707, 005 (2007). 
[23] I. S. Gradshteyn and I. M. Ryzhik, Table of Integrals, 

Series and Products (Academic Press, New York, 1980). 
[24] V. F. Mukhanov, JETP Lett. 41, 493 (1985); Sov. Phys. 

JETP 68 1297 (1988). 
[25] F. Finelli, J. Hamann, S. M. Leach and J. Lesgourgues, 

larXiv:0912.0522l [astro-ph.CO]. 
[26] G. Marozzi, Phys. Rev. D 76, 043504 (2007). 
[27] M. Lemoine, J. Martin and J. Yokoyama, Phys. Rev. D 

80, 123514 (2009); Europhys. Lett. 89, 29001 (2010). 
[28] D. Polarski and A. A. Starobinsky, Nucl. Phys. B 385, 

623 (1992). 



